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ABSTRACT 

Aims. In the framework of TRACE coronal observations, we compare the analysis and diagnostics of a loop after subtracting the background 
with two different and independent methods. 

Methods. The dataset includes sequences of images in the 171 A , 1 95 A filter bands of TRACE. One background subtraction method consists 
in taking as background values those obtained from interpolation between concentric strips around the analyzed loop. The o ther method is a 
pixel-to-pixel subtraction of the final image when the loop had completely faded out, already used bv lReale & Ciaravellal2006l 
Results. We compare the emission distributions along the loop obtained with the two methods and find that they are considerably different. We 
find differences as well in the related derive filter ratio and temperature profiles. In particular, the pixel-to-pixel subtraction leads to coherent 
diagnostics of a cooling loop. With the other subtraction the diagnostics are much less clear. 

Conclusions. The background subtraction is a delicate issue in the analysis of a loop. The pixel-to-pixel subtraction appears to be more reliable, 
but its application is not always possible. Subtraction from interpolation between surrounding regions can produce higher systematic errors, 
because of intersecting structures and of the large amount of subtracted emission in TRACE observations. 
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1. Introduction 



An important issue in the data analysis of coronal loops 
observations is background subtraction. Recent results ( 



The issue is critical especially in the analysis of obser- 
vations where the background is a significant fraction of the 
signal. This happens, for instance, in observation made with 



r-( ! |Del Zanna & Masonl l2003i iTesta et all 2002, Schmelz et al 



TRACE (e.glSchmelz et all I2003I IXschwanden & Nightingal 



2005I IReale & Ciaravella! I200oi Here we explore the depen 



2003, Aschwanden & Nightingale] I2005L IReale & Ciaravella 
2006& lAschwanden et aL I2008I) have established the impor- 



tance of separating the actual loop plasma from the diffuse fore- 
ground and background emission, that results from unresolved 
coronal structures and instrumental effects. The need of back- 
ground subtraction arises from the presence of many overlap- 
ping bright structures and of diffuse emission, n earby or along 



the li ne of sight, as well as from stray light ( DeFor est et al 



2009). The accurate extraction of the emission along the loop is 
necessary to apply standard diagnostic methods (like filters ra- 
tio) to derive physical quantities, such as temperature, or even 
to apply more detailed loop models. Wherever the background 
is non-negligible compared with the loop intensity, it would se- 
riously affect the extracted intensity along the loop. There is no 
standard and generally accepted method of background sub- 
traction; the procedure is then "operator- sensitive" and so are 
in turn the results. 



dence and sensitivity of the results on the background subtrac- 
tion method, by comparing two methods applied to the same 
loop, observed with TRACE in more than one filter band. We 
take advantage of the unique opportunity to apply two differ- 
ent and independent methods of backg round subtraction on th e 
same dataset. One of the methods (IReale & Cia ravella 2006 - 
hereafter RC06) has an objective support from the facts that we 
use as background an image of the loop region when the loop 
is absent, and that background variations are estimated to be 
small throughout the observation. 



In section 2 we describe the data analysis and background 
subtraction. Section 3 shows and compares the results with the 
two different background subtractions, including implications 
on temperature diagnostics with filter ratios. In section 4 we 
discuss the results. 
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Fig. 1. Loop region as observed with TRACE (256 x 256 pixels 
image) in the 171 A filter at the labelled time (lower panel). 
The loop under analysis is pointed with an arrow. The loop 
region is located in the inset of the Yohkoh/SXT full disk im- 
age (upper panel). The grey scale is inverted and linear for the 
TRACE image (< 8 DN/s/pix) and inverted and logarithmic for 
the Yohkoh image (between 10 and 150 DN/s/pix). 



2. Data analysis 

The Transition Region and Coronal Explorer (TRACE) is a 
NASA Small Explorer (SMEX) mission to image the solar 
corona and transition region at high angular and temporal reso- 
lution, operating since 1998 dHandv etalf l999). We analyze a 
TRACE observation of May 13 1998, with a 3.5htime sequence 
of 1024 x 1024 pixels full resolution image in two of TRACE 
filters (171 A , 195 A). This is the same data set selected and an- 
alyzed by RC06. The same 5 12 x5 12 pixels region of the whole 
field of view has been extracted for the analysis. From these 
data four images were selected, in each filter, at the following 
times: 06:36:57 UT, 06:59:35 UT, 07:39:26 UT, 08:29:34 UT 
for 171 A filter, and 06:37:18 UT, 07:00:06 UT, 07:39:47 UT, 



08:30:06 UT for 195 A filter. The filters have different sensi- 
tivity to temperature of observed plasma. Observations taken 
at the same time with different filters provide information on 
the plasma temperature. In particular, the ratio of the emis- 
sion in two different filter bands, compared with the ratio of 
the response functions of the filters, gives the temperature of 
the emitting plasma. Since the filter ratio may be altered by 
the diffuse emission, by the signals from other structures co- 
aligned along the line of sight and by the stray light, it is very 
important to estimate and remove such background emission. 
In TRACE observations the diffuse emission is typically very 
high and might represent most of the signal. 

2.1. Loop analysis 

We analyze the same loop as in RC06. The loop is selected on 
the TRACE images (FigJTJ, and appears as an entire loop in 
several 171 A filter images (Figj2j. It is clearly visible as well 
in the 195 A filter band (Fig|2]i. The loop is bright, i.e. ob- 
served with good count statistics and with a high contrast over 
the background. It was selected as far as possible free from 
other structures intersecting along the line of sight. The loop 
evolves in both filterbands. In the 171 A band it is faint ini- 
tially, then it brightens more reaching a peak of intensity around 
07:30 UT and then it fades again, and disappears at the end 
of the sequence (~ 10:00 UT, RC06). In the 195 A band the 
loop is brightest initially and progressively fades out. This loop 
evolution is even more clear in Fig. [3] which shows images of 
the loop region normalized to the final image at 10 UT, where 
the loop is not longer visible. Here we revisit the data at the 
four selected times. The data were treated with the standard 
procedures for TRACE data processing contained in the Solar 
Software (ssw) and the images were co-aligned with a standard 
cross-correlation technique. 



2.2. Background subtraction 

Here we compare the method already adopted in RC06 
with an alternative and indepedent approach, based on in- 
terpolation between emission values in a region close to 
(but outside) the loop, sim i lar to t hose in Testa et al. 2002, 
Aschwanden & Nightingale] l2005i Schmelz et alj 2003 and 
Aschwanden et alJl2008l As a first step, we measure the emis- 



sion of the loop. We define a strip enclosing the loop in both 
TRACE filter passband, down to the visible footpoints and even 
beyond them, and divide it into sectors, as shown in Fig]4] We 
have analyzed strips of different widths. A width of 10 pixels is 
a good compromise between too low statistics and overcoming 
too much the loop borders. With this choice we end up with 
27 similar almost square sectors (RC06). Then, we mark two 
other strips, parallel and concentric to the one used to extract 
the loop emission, one outer and one inner of the imaginary 
circle which the loop is part of. We divide these two additional 
strips into sectors aligned to those of the central strip (Fig J4j. 
These two other strips are used to extract the emission for back- 
ground subtraction. We put them as close as possible to the 
loop of interest, avoiding other structures near, although dis- 
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Fig. 2. Loop region (top row: 171 A ; bottom row: 195 A ) at subsequent times (from left to right): 06:36:57.000 UT, 
06:59:35.000 UT, 07:39:26.000 UT, 08:29:34.000 UT for 171 A and 06:37:18.000 UT, 07:00:06.000 UT, 07:39:47.000 UT, 
08:30:06.000 UT for 195 A. The grey scale is inverted and linear between 0.4 and 7 DN s pix 1 for all images. The frame in 
the top right image is the part shown in Fig. [3] 



tinct from, the analyzed loop. Both strips also have a width of 
10 pixels, and are divided into 27 sectors, so to have a one-to- 
one correspondence of these sectors with the sectors patching 
the loop. In each of the 27 sectors of the two external strips we 
compute the mean emission value per pixel. To each sector of 
the central strip we assign as background the emission value 
linearly interpolated between the two corresponding sectors of 
the external strips (Fig|4ji. 

For the sake of completeness, as a further refinement we 
have considered also a least-squares quadratic interpolation. 
We have then considered two further concentric strips around 
the analyzed loop, one outer and the other inner. However, the 
additional outer strips are so far from the central loop, that they 
end up to include other bright structures which severely con- 
taminate the procedure and locally alter the otherwise smooth 
emission distribution. Therefore, in our case the higher order 
interpolation does not lead to an improvement of the results. 
Anyhow, we point out that our basic linear interpolation is 
made between values averaged over a sector, and therefore is 
already and intrinsically more accurate than a linear interpola- 
tion made between values in single pixels. 

We remind that the method of background subtraction of 
RC06 is based on the fact that the loop disappears at the end of 
the image sequence. The last image (around 10 UT) is then sub- 
tracted pixel-by-pixel from all other images, assuming that the 
structures surrounding and crossing the loop of interest along 
the line of sight do not change much during the observation se- 
quence. As an estimate of the fluctuations of the background 



during the observation, RC06 measured an average pixel-by- 
pixel standard deviation of 13% both in the 171 A and in the 
195 A filter band. 

A more accurate subtraction would be obtained if the sub- 
tracted image resulted from the interpolation between an image 
before the loop switches on and one after the loop switches off. 
Unfortunately, this loop has a very long lifetime and in order to 
find images before the loop ignition we have to go back more 
than 25 ks earlier than the start time of our analysis. In our opin- 
ion, this time is too long to pretend that the environment did not 
change in the meantime and to have a reliable background, and 
we prefer to keep the final image as the only image for subtrac- 
tion. 

We also point out that the image to subtract was taken at 
10:00 UT, significantly apart in time from the analyzed images 
taken between 06:30 and 08:30 UT. The perspective (i.e. view- 
ing angle) of the observed loop was changing during that in- 
terval and this may impact the RC06 scheme which assumes 
that each line of sight intersects the same column of plasma 
throughout the entire sequence. However, we do not expect a 
significant effect since we estimate a change of the viewing 
angle of the order of 1 - 2° or less, as confirmed by a visual 
inspection of the images after subtraction. 

Thes e are not the only t wo methods of backgr ound sub 
traction. 



Schmelz et ai]l2007l andfSchmerz et al. 2003 used two 



other lower-order methods of background subtraction: ( 1 ) con- 
stant background subtraction; (2) pixel pair background sub- 
traction. In the former method they subtract a constant value 
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Fig.3. Zooms of the loop region in Fig. [2](frame in the top right image of Fig. |2j. The images (top row: 171 A ; bottom row: 
195 A ) are normalized to the final image of the sequence (at 10:00 UT) where the loop is no longer visible, to enhance the loop 
evolution. The grey scale is linear between 0.4 and 2.5 (1 means that the emission value is the same as that of the final image). 



extracted from a single background pixel chosen above the loop 
apex; in the latter they selected a background pixel for each 
loop pixel. 

3. Results 

3.1. Loop Emission 

Fig [5] and [6] show examples of the application of our interpola- 
tion method of background subtraction. The lines join the data 
points, obtained as mean values of the pixel emission in a sec- 
tor. The error of each mean value is estimated as the standard 
deviation of the mean in that sector. The error has been conser- 
vatively rounded to 0.05. The lower panels of Fig(5]and|6]show 
the results of the subtraction of values interpolating between 
the dotted and the dashed curves, from the solid curve. 

The loop emission is significantly lowered by the back- 
ground subtraction, in some points even completely cancelled, 
due to the overlap with other bright structure close to, but dis- 
antangled from our loop. For a direct comparison of this and the 
RC06 background subtraction methods, we now plot the emis- 
sion along the loop obtained with the two methods in the same 
figure. Fig|7] shows the results in the 171 A filter and in the 
195 A filter. We can see that the emission is significantly dif- 
ferent after the background subtraction with the two methods: 
the profiles mostly differ in the central part, which is predicted 
to be very faint at any time after background subtraction by 
interpolation method. On the contrary, with the latter method 
the footpoints are found to be mostly brighter than with the 
pixel-to-pixel subtraction, although we see some agreement at 



some times (e.g. right end at 06:36 UT and 07:00 UT in both 
filters). While the result regarding the footpoints may be de- 
bated, since the footpoint regions cannot be clearly resolved 
in the original data (Figf2]i, there is no doubt that the bulk of 
the loop is clearly visible at many times (Figf2] and FigO, 
and that, therefore, the emission cannot be so low, in contrast 
with interpolation method and in agreement with pixel-to-pixel 
method. Therefore, in the case where the two methods most 
differ, RC06's pixel-to-pixel method surely provides a more re- 
liable result. There is some agreement between the profiles at 
the latest time. 

Moreover, we notice in Fig|7] that, since the subtracted 
emission obtained with interpolation method is invariably very 
low except near the loop footpoints, it does not show a clear 
evolution. Instead RC06 obtain that the background-subtracted 
central section of the loop has an evolution in agreement with 
the evolution observed in the images (Figf2]and Figf3]i. Again, 
therefore the results obtained with the method of RC06 appear 
more reliable. 

3.2. Temperature diagnostics 

To further assess the differences between the background sub- 
traction methods, we have explored the temperature diagnostics 
obtained with filter ratio. It is well-known that for an optically 
thin plasma, isothermal along the line of sight, the ratio of the 
emission detected in two different filters is a function of the 



temperature only (e.g. lRosner et al.l 19781) . From the ratio value 
measured in an image pixel we can then derive a temperature 
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Fig. 4. Background subtraction with interpolation method. 
Upper panel: the strips for the loop analysis (central) and for 
the background subtraction (outer and inner) are marked in the 
loop region (07:39:26.000 UT, 171 A filter). Lower panel: 
Emission profile (thin solid line) measured on a cut across the 
loop (marked in the upper panel). We compute the average 
emission (thick horizontal lines) in the corresponding sectors 
intersected by the outer and inner strips (bounded by vertical 
dashed lines). The background value (central horizontal line) is 
computed by linear interpolation of the average emission val- 
ues of the outer strips. This is then subtracted from the total 
emission measured in the central strip (upper solid line). 



value in that pixel, and build therefore a proper thermal map. 
In the case of the narrow band 171 A and 195 A filters of 
TRACE, the relationship between temperature and filter ratio 
195/171 is monotonic onl y in a certain temperature range, i.e. 



0.7 < T < 1.8 MK (e.g. lAschwanden et al.ll2000l) . which is 
also the range where the filters responses are the highest. So it 
is reasonable to assume that the detected plasma is all in that 
temperature range and this allows us to measure the tempera- 



ture along the loop, after background subtraction. As an exam- 
ple, Fig|8]shows the filter ratio 195/171 and corresponding tem- 
perature along the loop at 07:40 UT computed after the back- 
ground subtraction with the pixel-to-pixel method (also shown 
in RC06) and with the interpolation method. It is immediately 
apparent that the ratio profile obtained with the latter method is 
much more irregular, so that we cannot define an overall trend. 
The other method instead leads to a more regular profile, co- 
herent with the clear visibility of the loop in the filter ratio map 
(RC06), in spite of the dip on the left leg due to another cross- 
ing bright structure. We find similar trends in the corresponding 
temperature profiles along the loop. 

We have to point out that, although both profiles indicate an 
overall small temperature variation along the loop, the pixel- 
to-pixel background makes the difference in identifying vs not 
identifying a trend. The globally fla t temperature pr ofile is a 
typical result from TRACE data (e.g. lLenz et al . 1999) and has 



been explained partly with an instrumental bias (IWeber et al 
2005). However, RC06 pointed out that the filter ratio diagnos- 
tics can be anyhow meaningful and is sensitive enough to de- 
tect the loop progressive cooling. Fig.|9]shows the loop average 
filter ratio obtained after the two different background subtrac- 
tions at four different times. From the evolution obtained with 
the interpolation method it is much less clear that the loop is 
cooling, confirming that the sensitivity of the temperature di- 
agnostics is higher with the pixel-to-pixel background subtrac- 
tion. 

4. Discussion 

The background subtraction is important for the analysis of 
TRACE coronal loops: the background signal is very high 
and influenced by the presence of many bright structures near 
and perhaps entangled with the analyzed loop, along the line 
of sight. Recent resul ts ( Del Zanna & Mas onll2003llTesta et al.l 
| 2002[ ISchmelz et al]|2003l lAschwanden & Nightingale! 120051 
Real e & Ciarav ella 2006, Aschwande n et al.1 2008) have estab- 
lished the importance of separating the actual loop plasma 
from the diffuse foreground and background emission that re- 
sults from unresolved coronal structures and instrumental ef- 
fects (stray light). Here we catch the unique opportunity to 
compare two different and independent methods of background 
subtraction: one that subtracts the interpolated emission of two 
off-loop strips; and the other that subtracts pixel-by-pixel the 
complete image after the disappearance of the loop (RC06). 
Extracting the background emission with two strips as near as 
possible to the loop of interest, we certainly consider and re- 
move emission from structures that can intersect transversally 
the loop. This is true as long as the "contaminating" loops run 
across the target loop, while it is not the case if this happens at 
an oblique angle. 

We have shown that not only do different methods lead 
to different emission profiles, but also the subsequent diag- 
nostics is affected. With the pixel-by-pixel subtraction we are 
able to derive quite coherent filter ratio profiles along the 
loop and quite coherent temperature evolution: the loop is 
globally cooling as expected from the sequence of appear- 
ance/disappearance in the different TRACE filters. Such infor- 
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Fig. 5. Emission along the loop in the 171 A filter, at 06:59:35.000 UT (left column), and 07:39:26.000 UT (right column). Upper 
panels: emission along the loop strip (solid line), outer strip (dotted line), and inner strip (dashed line). Lower panels: emission 
after background subtraction with the interpolation method. 
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Fig. 6. As Fig|5]for the 195 A filter, at 07:00:06.000 UT (left column), and 07:39:47.000 UT (right column). 
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Fig. 7. Emission along the loop at the labelled times in the 171 A filter (left column) and in the 195 A filter (right column). 
All panels show the emission along the loop, after background subtraction with two different methods: pixel-by-pixel method 
(RC06, dashed line), and interpolation method (solid line). 



mation is much less clear after using the other method of back- 
ground subtraction. 

The RC06 background subtraction looks therefore a more 
reliable method. It is in principle very accurate, if the loop of in- 
terest is the most variable structure in the field of observation; it 
is direct, at variance from the interpolation method; it is applied 
pixel-by-pixel, allowing us to derive "background-subtracted 
images" and therefore to have a visual feedback, and to ana- 



lyze all loop pixels, instead of sampling them at selected posi- 
tions. On the other hand, this method works well as long as the 
structures surrounding and crossing the loop of interest along 
the line of sight do not change much during the observation. In 
other words, the RC06 method cannot take time-variations of 
the background emission into account and we cannot exclude 
that crossing structures vary during the observation. We also re- 
mark that the pixel-by-pixel subtraction could be applied only 
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Fig. 8. Filter ratio 195/171 (top) and corresponding tempera- 
ture (bottom) along the loop at 07:40 UT computed after the 
different background subtractions, pixel-to-pixel (dashed) and 
interpolation (solid). The arrow is an upper limit among inter- 
polation data. 



because the loop disappears at the end of the image sequence. 
This condition can be matched only by evolving loops. For 
loops keeping steady during the entire observation, other meth- 
ods must be use d. The interpolation of out-loop emission has 
been often used (Testa et al . 2002, Aschwanden & Nightingale 



20051 ISchmelz et al J I20Q3L and lAschwanden et all 12008). but 
according to our analysis, it may lead to severe systematic er- 
rors because it may be affected by other structures close but 
distinct from the loop under analysis. 

Fig- HI points to the difficulties encountered by the interpo- 
lation method since it shows that relatively close to the target 
loop there is another structure with similar intensity. Although 
the interpolation method should provide a good estimate of the 
background around the two loops, it is not able to take the over- 
lap of the two loops into full and proper account. A scheme 
which fits the two loops separately together with a background 
estimated as above may represent a better solution. 

Other methods of background subtraction are instead too 
"operator-sensitive", because they are based upon the meticu- 
lous selection of single background pixels, intentionally avoid- 



Fig. 9. Evolution of the loop average 195/171 filter ratio, sam- 
pled at the four times analyzed in this work and obtained after 
the different background subtractions, pixel-to-pixel (stars) and 
interpolation (triangles). 

ing in this way structures that can cross the loop of interest and 
can alter significantly the loop emission. 

In conclusion, this work confirms and qualifies how a re- 
liable background subtraction is a delicate and difficult task 
for the analysis of coronal loops observed with TRACE. The 
problem can be of course greatly reduced in observations with 
different instruments where the amount of instrumental back- 
ground emission is lower and where therefore the presence of 
systematic effects has less influence. 
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